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Abstract: We consider particles in three-dimensional space, which have a certain 
probability to find themselves in a thin layer ( "plane" ) , where they are assumed to be well 
described by a planar Hamiltonian and are subject to Aharonov-Bohm-type interaction. 
We demonstrate that their planar motion is then anyonlike, with the "effective statisti- 
cal parameter" proportional essentially to the square of the probability. We also show 
that charge-flux composites with solenoids of finite length, provided they are themselves 
fermions, can form a bound state in which they behave like anyons, without any external 
potential confining them to a plane. 
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1. As is well-known, fractional statistics is possible only in two-dimensional space |l]] 
(in one-dimensional space as well, but this is not our concern here); this is dictated by 
general topological arguments. On the other hand, the real world is three-dimensional, and 
two-dimensionality may only be regarded as an approximation. Therefore for fundamental 
particles, fractional statistics is impossible. On the other hand, it is intuitively clear that 
if a particle is confined in one direction strongly enough, a two-dimensional model should 
describe its behavior well enough. Thus, a question arises about a quantitative description 
of the transition from a "purely three-dimensional" regime, with no anyons possible, to a 
"quasi-two-dimensional" one, where they may be considered. 

Since no particle itself can be an anyon, in reality it is only possible to obtain fractional 
statistics effectively by means of some interaction. The most relevant example seems to be 
the model [§] where the "particles" (more generally, quasiparticle excitations) are made 
of charges and infinitely thin solenoids and therefore interact a la Aharonov-Bohm. If 
the solenoid has an infinite length, then nothing depends on the third coordinate (z) and 
the two-dimensional description is exact. If it is much longer than the domain of motion 
along the z axis, then this description may be regarded as a good approximation. We are 
going to clarify what happens if one goes beyond that approximation. We will consider 
the two-particle problem, including the calculation of the second virial coefficient, and 
the problem of the A^-particle ground state. In the second part we will demonstrate that 
under certain conditions, anyonic interaction itself may confine the particles to a plane, 
forming a bound state. 

Anyons are particles that live on a plane and whose two-particle Hamiltonian is 

H(R,<S>;r,<p) = H CU (R, + H (r, <p) + AH (r, cp; a) , (1) 
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where Hqm(R,$) is the (irrelevant) center-of-mass Hamiltonian and 

• . Id 2 Id Id 2 muj 2 9 . , 

H (r, if = —J ~ J + ' 2 

m or A mr or mr z oip z 4 

. rr , . 2ia da 2 

AH (r,<p = —jt + —- 3 
mr z ocp mr z 

Here the so-called regular gauge [|| has been chosen, in which the wave function is sym- 
metric (or antisymmetric) with respect to interchange: \P(r, — (p) = ^(r,(p) [— ^f(r, (p), 
respectively]; a is the statistical parameter, and a harmonic attraction has been added by 
hand in order to discretize the spectrum. The nature of the interaction (^), for example 
in the context of the fractional quantum Hall effect, may be understood as follows: The 
particles are confined to a plane (the interface between two semiconductors etc.) by an 
external potential and being on this plane, they effectively have a magnetic flux "glued" 
to them; it is then the charge-flux interaction that is described by AH(r, <f>) (it changes 
the kinetic angular momentum by a). Imagine now a situation in which the "plane" has 
a finite thickness^ 21 while the particles can move in the z direction within a range of 
width 2D; when a particle is within the "plane", it acquires the flux (stretched in the z 
direction) and therefore interacts with the other particles in the above-described manner, 
otherwise it does not interact with them. Clearly, the relevant two-particle Hamiltonian 
takes the form 

H = H ± (z 1 ) + H ± (z 2 ) + H CM (R, *) + H (r, <p) + 9(1 - \ Zl \)8(l - \z 2 \)AH(r, <p; a), (4) 
where 



lm dz' 

V(z) is some confining potential, which determines D, and 6(z) is the step function. 



3 It is to remind of this thickness that we use the quotation marks. 
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This model can be said to describe an effective change of statistics in three dimensions. 
Indeed, if D < I (to be more exact, if V(z) = oo for \z\ > I), we arrive to the above- 
discussed situation of infinite solenoids, and the planar motion decouples from the z motion 
and is purely anyonic. On the contrary, if D S> /, the last term in (Q) has essentially no 
effect, and one has pure bosons (or fermions). In the most interesting intermediate case 
D > I, the problem does not allow for a separation of variables and therefore in general 
is not solvable. However, if \a\ <C 1, then it is possible to consider the last term as a 
perturbation. The unperturbed solutions are (here and after we completely ignore the 
planar center-of-mass excitations) 

^Sfeain = V^l)V>fc 2 ( z 2)^Ln(r, (p), (6) 

where ip^r(z) are the eigenfunctions of H and ipLnif,^) are the eigenfunctions of Hq; in 
particular, for Hq as in (|2|) we have 

Mr, <p) = Cfor^ (-n, \L\ + 1; -^r 2 ) exp (-^r 2 ) , (7) 

and correspondingly 

41l„ = < + Ei + (2n + \L\ + l)u (8) 

(L is even/odd for bosons/fermions; in what follows we will consider bosons, until otherwise 
stated). Now, for the two-anyon Hamiltonian, first-order perturbation theory gives [Q, ||] 
the result which is in fact correct for any a (l|, U, 

au), L > 0, 

{^Ln(r,<p)\ AH(r,ip;a)\ip Ln (r,(p)) = { \ a \u, L = 0, (9) 

-aw, L < 0. 

Hence the first-order energy in our problem is easily found to be 

4!Ln = E i + E k + ( 2n + \ L + w kl w k2 a\ + (10) 
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ti(z) dz (11) 



is the probability that a particle with the wave function i(^k( z ) stays within the "plane". 

This result is natural enough, illustrating that the particles are anyons "inasmuch 
as they are in the plane": The planar part of the energy corresponds to that of two 
anyons with an "effective statistical parameter" equal to a times the probability that 
both particles are in the "plane" . Passing from D <C I to D ^> I will correspond to passing 
from Wfa 2 ~ 1 to Wfa 2 ~ 0. Note by the way that for perturbation theory to be applicable, 
it is in fact not necessary that \a\ <C 1, it is sufficient that Wk 1 Wk 2 \a\ <C 1 [the result (||) 
is exact anyway]. 

As an illustration, let us calculate explicitly the "planar" second virial coefficient in 
this situation, defined as usually 0, 



(f) 2 Ht)< 
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A = y/27r/mT being the thermal wavelength. (Here and further, Z stands for three- 
dimensional partition functions, Z and Z 1 - for planar and perpendicular ones, respec- 
tively.) Let the perpendicular potential be an infinitely deep well: 

, 0, 1*1 < D, 
V(z) = (13) 



Then 



where 



oo, \z\ > D. 
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One has 
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where 
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and 9i(z, q), i = 1, . . . ,4 are elliptic theta functions. On the other hand, 



(17) 
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Zi (19) 



(for L = 0, the last term in the square brackets has to be replaced by Wk 1 Wk 2 \ a W)-, an d 
expanding to the first order in a, one eventually gets 

1 



+ H 2 \a\ A 2 , 



(20) 
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^ oo 

» = ^xE^ ex p(-^/ T ) 



y l fc=l 



7 



1 + 



1 - j*o 6» 4 (vr7x/2, g) cte 



(21) 



3 (0,<?) -1 

That is, /i is the thermal average of the quantity w^, and again the "effective statistical 
parameter" is proportional to fi 2 , reflecting the fact that two-particle interaction is in- 
volved. Clearly, for q = (which corresponds to T = 0) one has n = u>\ = 7 + sin(7T7)/7r, 
while for q = 1 (T = 00) the probability distribution along z becomes uniform and (i = 7. 
Fig. 1 shows the dependences ^(7) for these two limit cases (dashed lines) as well as for 
an intermediate case q = 0.75 (solid line). 
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Let us now extend the above considerations to the iV-particle case. The complete 
iV-anyon spectrum is not known, for N > 3, but the ground state for small enough a is 
known: namely, if 

"'M-S^ < 22 > 

is the one-particle (planar) Hamiltonian, H n (ri, . . . , r n ; a) is the n-anyon Hamiltonian in 
the regular gauge, and 

n 

AH n (r 1 , . . . ,r n ;a) = H n (r 1: . . . ,r n ;a) -J^H^rj), (23) 

i=i 

then the first-order correction from AH n to the ground state of J2j H\ (r j ) , with the wave 
function Vo(ri, • • • , r n ) = c n exp YTj=i r j) , is (see g again) 

Ttifl — 1 ) 

(V>o(ri,...,r n )|Afr n (ri,...,r n ;a)|^o(ri,...,r n )) = \a\u. (24) 

What will be the generalization of Eq. (Q) for this case? Introduce a quantity 

= °( l - M) '"0(1- \z k J)0(\z kn+1 \ -/)••• 9(\z kN \ - l), (25) 

where k\, . . . , k n are any n different numbers from the set {1, . . . , N} and k n+ \, . . . , kjy are 
the remaining ones from this set. Thus, kn is equal to 1 if and only if the particles with 
numbers ki, . . . , k n are in the "plane" and all the others are out, otherwise it is equal to 
0. Following our assumption that the particles interact only when they are in the "plane", 
we get 

N N N 

H N = Y,H ± (z J )+J2Hi(r J ) + J2 E ^L.k n AH n(rk 1 ,---,r kn ;a), (26) 
J=l j=l n=0 {ki...k n } 

where the last sum is performed over all possible choices of n numbers k\, . . . , k n , for a 
given n; the number of terms in this sum is the binomial coefficient Cjy. Assuming again 
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that it is possible to consider the last term as a perturbation, the first-order ground state 
energy of ( p6|) will be 

n r _ -n 

E$ = NEq +Nuj + Y: C n N w n (l - Wo ) N -"^— L l\ a \ u 

n=0 1 

= NE^ + Nto+ N ( N ^~ ^ w g| a |cj, (27) 

where again wq is the probability for a particle to be in the "plane", as in (p]). Once 
again, this corresponds, with the subtraction of the z motion, to the ground state energy of 
N anyons with an "effective statistical parameter" WqCx; this is again natural, taking into 
account that upon introducing the ansatz necessary to make perturbation theory work || , 
the interaction Hamiltonian becomes a sum of pair terms. Alternatively, if N 3> 1, one 
may represent the last term as ( w °^ \a\u, interpreting it as the ground state energy of 
woN anyons, woN being the average quantity of them in the "plane" . 

2. So long we assumed the particles to be confined to a fixed "plane" by an external 
potential, as it happens, for example, in the quantum Hall effect; then the whole spectrum 
is certainly discrete. Now we are going to show that a two-particle bound state can as 
well emerge even when the confining potential is absent, so that the z motion is free; in 
this situation there is no fixed "plane", and we suppose that the interaction is present 
for \z\ — 22 1 < I. This would correspond to the particles being charge-flux composites 
with solenoid length I and the function <£(r, z), the flux felt by one charge from the other 
solenoid when their centers are separated by the vector (r, z), being approximated by 
(j> ■ 9(\l\ — z); dropping the r dependence is justified at least for r <C I, when the flux lines 
close far away from the charges. Our idea is simple enough: If the particles themselves are 
fermions, then the energy of the two-particle planar ground state is 2uj outside the "plane" 
but (2 — \a\)co inside (note that here and further, a will mean what (1 — a) usually means, 
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i.e. the deviation from Fermi statistics); now, in three dimensions a state with energy E 
such that (2 — \a\)u < E < 2uj can exist — the z dependence of its wave function will be 
a superposition of plane waves inside the "plane" , but exponential damping outside, due 
to the right inequality. (See, e.g., || for a discussion of an analogous situation, when a 
bound state exists in a classically unbound system.) 

For an explicit computation, recall the general case: Let a system have a boundary 
z = I and thus be characterized by a Hamiltonian of the form 

[ H-(r)+H-(z), z < I, 

H(r,z) = { K1 {h (28) 

I H+(r) +H+(z), z>l. 

Suppose that the following information is available: (i) the solution of the spectral problem 
for H~(r) and H + (r), that is the energies and wave functions in 

F ± (r)^(r)=£:±0±(r), (29) 

and (ii) the wave functions ip {£, z), for every £, such that 

H i: {z)il) ± {S, z) = £tp ± {£ 1 z) (30) 

and that Jq ±0 ° 1^(8, z)\ 2 dz is finite. [If £ takes a value within the discrete set of eigen- 
values of the extension of H + {z) to z E (— oo, oo), then the integral /f^ \ip + (£, z)\ 2 dz is 
finite as well, for any other £ it is of course divergent, the same for ip~(£, z).\ Then it is 



possible to solve the problem (28). Indeed, since in each of the two regions the variables 



separate, the wave function has to be searched for in the form 

*(r,z) = { " (31) 
Ec+^+^-f+z), z>l; 
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by construction, it satisfies the equation H(r, z)^(r, z) = E^(r,z) and is normalizable. 
The two boundary conditions take the form 



(32) 



and the same with V ;± replaced by if)f = difj^/dz. Scalar ly multiplying both these condi- 
tions by (ffc (r), we get 



,0 = E c+$ fcm v 2 + (£ - s+, i), 



where 



<f km = J fc -(r)0+(r)d 2 
is the overlap integral, and finally, excluding c^T , 



(33) 



(34) 



(35) 



where 



ijA(E - £+, l)^(E - £-,l) - tf(E - £+, l)r (E - £~,l) 



km- 



(36) 



The equation det A^n = then determines the eigenvalues of E. 

Let us calculate the ground state in the problem at hand. Symmetry with respect to 
z = is present, hence the ground state wave function has to be even with respect to z 
(nodeless), which means it is enough to consider z G [0, oo); we have 



H~{z) = H + (z) = 0, 

H+(r) = H (r,<p), 

H-(r) = H (r,<p) + AH(r,<p;a); 



(37) 
(38) 
(39) 
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in view of the aforesaid, one has to choose 

ip + (£,z) = exp(-V-2m£z) (£ < 0), (40) 

^~(£,z) = cos(v / 2^£z) {£ > 0). (41) 

Again, let us restrict ourselves to the simplest case \a\ <C 1, whence (&km — 5km + ctXk m \ 
since all terms in det^4/t m , except rifc^fcfc) will then have at least a factor a 2 , we can 
neglect them. We are to search for E in the form 

E = (2 - k\o\)lo, (42) 

with < k < 1, the only possibility is then to have A 00 = 0; now £q = (2— \a>\)u, £q = 2u, 
and an equation for k follows 

tan yj2\a\£(l-K) = <f^, (43) 

where 

i = mtul 2 . (44) 

The asymptotic expressions are k ~ 2\a\£ for |a|^ <C 1 and k —>■ 1 for £ — ► oo; the whole 
dependence «(£), for \a\ = 0.5, is sampled on Fig. 2. Thus, in this situation it is the ratio 
of the solenoid length (I) and the planar scale (l/^mw) that matters: The first one being 
much greater brings us again to the case of ideal anyons (E — > fj~), in the second one, 
clearly, the length of the "tail" of the wave function (along z) is much greater than I, and 
E — > <S 1 ', which means that the particles become essentially free; however, the bound state 
exists always, despite the absence of a confining potential. (Note that the above-discussed 
replacement of <J>(r, z) by the step function would in fact be good in the first regime only, 
i.e. for £ 3> 1; however, the qualitative picture will remain the same as described in any 
case.) 
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We come to the conclusions that: 

1. If the particles (generally speaking, not pointlike) have a finite probability to find 
themselves within a "plane" — a region where they are subject to anyonic interaction, they 
may be regarded, in what concerns their planar motion, as anyons with the statistical 
parameter proportional typically to the square of this probability This is a natural way 
of realizing how one can (effectively) get fractional statistics within an underlying three- 
dimensional model; 

2. For charge-flux composites with solenoids of finite length which are fermoions by 
themselves, the anyonic interaction, having a character of attraction, can lead to formation 
of a bound state even if their longitudinal motion is by no means restricted. 

We wonder whether a specific microscopic model could be found in which this descrip- 
tion would work properly In particular, it would be an interesting problem to study in 
more detail the influence of possible "jumps" of particles (excitations) off the plane, for 
example, on the fractional quantum Hall effect, in particular to observe how the effect 
disappears when confinement to the plane becomes less strong. 
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FIGURE CAPTIONS 

Fig. 1. The probability ^(7) that a particle stays within the "plane", for three cases: 
q = (dashed curve); q = 0.75 (solid curve); q = 1 (dashed straight line). 

Fig. 2. The function k(£) which determines the energy of the bound state, for \a\ = 0.5. 
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